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We develop a theory for the electron-phonon interaction effects on the electronic properties of 
graphene. We analytically calculate the electron self-energy, spectral function and band veloc- 
ity renormalization due to phonon-mediated electron-electron interaction. We find that phonon- 
mediated electron-electron coupling has a large effect on the graphene band structure renormaliza- 
tion, and our analytic theory successfully captures the essential features of the observed graphene 
electron spectra in the ARPES experiments, predicting a kink at ~ 200meV below the Fermi level 
and a reduction of the band velocity by ~ f — 20% at the experimental doping level. 
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Electronic properties of graphene are of fundamental 
importance because of its linear two-dimensional Dirac- 
like energy dispersion, attracting considerable current at- 
tention [ij. Although there has been a great deal of 
theoretical work on the effects of disorder and electron- 
electron (e-e) interaction on graphene electronic prop- 
erties, the effect of electron-phonon (e-ph) interaction in 
graphene has not been studied extensively. In this Letter, 
we present a leading order many-body theoretic analysis 
of e-ph interaction-induced renormalization of graphene 
electronic properties, comparing our theoretical results 
critically with existing experimental results . The mo- 
tivation of our theoretical work comes from the beautiful 
recent experimental ARPES (angle-resolved photoemis- 
sion) study 0] of graphene electronic spectrum, and our 
theory explains the main features of the observed ARPES 
spectra as arising from the e-ph interaction. 



We focus our investigation on the effect of e-ph inter- 
actions on band structure renormalization by calculating 
the electron self-energy due to the phonon-mediated e-e 
interaction. Graphene has a two-dimensional honeycomb 
real-space lattice structure comprising two interpenetrat- 
ing triangular sublattices A and B, which translates to 
a reciprocal-space honeycomb structure with the hexag- 
onal Brillouin zone cornering at the high-symmetry K 
points. In the vicinity of these K points (so-called Dirac 
points), the low-energy excitations have a linear energy 
spectrum described by the low-energy bare Hamiltonian 
Hq = vcr ■ k, where v w lO^ms"^ is the Fermi veloc- 
ity of the Dirac fermions, <t is the set of Pauli matrices 
representing the two (A and B) sublattice "pseudo-spin" 
degrees of freedom (we let ft = 1 throughout unless it 
is written out explicitly). This Hamiltonian describes a 
cone-like linear energy spectrum with ekx — Ae^, where 
ek = vk, X is the chirality label standing for the con- 
duction band (A — 1) or the valence band (A = —1). 
The e-ph interaction was originally derived in Ref. [sj, 
from which we can write the second-quantized form in 



the momentum space as: 
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where Ck = [ak bk]"^ is the two-component electron anni- 
hilation operator in the momentum space for the A and 
B sublattices and dq is the phonon annihilation opera- 
tor. Phonon modes in graphene couple neighbouring A- 
sublattice and B-sublattice carbon atoms through bond 
stretching and bending, so that the e-ph coupling be- 
comes an off-diagonal matrix in the pseudo-spin space 
[i] gM, where the constant g = - {/3v /b'^ ) ^h/2NMcUJo 
gives the magnitude of the e-ph coupling, with loq = 
0.196eV the optical phonon frequency for graphene from 
Raman scattering experiments [H, IJ, N the number of 
unit cells, Mc the mass of a carbon atom, b = a/y/3 
the equilibrium bond length between adjacent carbon 
atoms and /3 — dln7o/dlnfe 2 is a dimensionless pa- 
rameter that gives the change of the nearest-neighbour 
tight-binding matrix element 70 with respect to the bond 
length b [||. The matrix M, for LO/TO phonons, is 
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with Mab = — 1 or i and M-qa = 1 or z for LO or TO 
phonons, respectively, and (j>q = tan^^ {qy/qx) the azith- 
muthal angle of the momentum q. We use this model 
e-ph interaction in our analysis. 

In the leading order theory the e-ph coupling Vcp could 
be transformed through a canonical transformation to an 
equivalent phonon-mediated e-e interaction V^J^. It fol- 
lows from Eqs. ((ij-llll) that the e-e scattering via phonon 
emission/absorption will flip the sublattice label from A 
to B or vice versa, and we obtain the following expression 
for the phonon-mediated e-e interaction resulting from 
the e-ph interaction of Eq. ([T]) (we adopt the convention 
where summation over repeated indices is implied): 

V,f = y' 'D\q,T2-Ti)M^,^{q)Mp,p{-q) 
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'cl^-ql3'CkiaCk2f3, (3) 
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sum in Eq. ([6]) can be performed explicitly at zero tem- 
perature. For concreteness, we take the case of n-doped 
graphene and evaluate the self-energy for quasipar- 
ticles in the conduction band (for p-doped material, one 
should consider Efc_ for quasi- holes in the valence band). 
The quasiparticle decay rate is proportional to the imag- 
inary part of the self-energy, which can be evaluated as: 
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FIG. 1: (Color Online) The real and imaginary parts of the 
self-energy at k — kp for density n = 10^"^ cm^^. The red 
dashed line shows the position of the Dirac point. 



where V^^q^Ti — T2) is the non-interacting retarded 
phonon Green function in the time domain. In the fre- 
quency domain, it is given in terms of the Matsubara 
bosonic frequency iqn as 'D'^{q,iq„) = 2uJo/[{iqn)'^ — loq]. 
Eq. ([3]) is the central quantity of this paper from which 
other quantities such as self-energy are derived. To find 
the self-energy for electrons in the conduction band/holes 
in the valence band, we work in the chiral basis (i.e. the 
diagonal basis of the Hamiltonian Hq) and denote quan- 
tities in the chiral basis with an overhead tilde. We shall 
focus ourselves only on LO phonons, as the calculations 
for TO phonons parallel that for LO phonons. The ma- 
trix elements for electron scattering from chirality A — > A' 
and fi fi' through phonon emission/absorption are 

{k + qX'\M{q)\kX) = i[A'e-*('^''+«-'^'') - Ae'^'^"""*')], 

{k-q^l'\M{-q)\k^i) = i[-^'e-'('^^-'-'^') +^e*(^''-^')], 

(4) 

and the corresponding phonon-mediated e-e interaction 
in the chiral basis can be written as 

Vt^lg' E T^°iq,r,-r2){k,+qX'\M{q)\k,\) 
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where the operators c\ c in Eq. ([5]) denotes the creation 
and annihilation operators in the chiral basis, with X, fi = 
±1 the chirality. 

The effective many-body velocity/mass renormaliza- 
tion to the band structure comes in large part from the 
electron-phonon interaction with the Coulomb interac- 
tion yielding a quantitatively small correction. The ef- 
fective velocity renormalization due to screened Coulomb 
interaction was considered in Ref. [^. In this Letter, we 
shall focus on the many-body effects of the e-ph inter- 
action without the effects of Coulomb interaction. From 
Eqs. (lll-lO, we have derived the following expression 
for the electron self-energy in the chiral basis due to the 
phonon-mediated e-e interaction Eq. ([5|): 
2 

X q,iqn 

[1 T Acos(0fe+g - 2(l)g)] , (6) 

where G^)^{ikn) — l/{ikn — £,k\) is the non- interacting 
electron Green function in the chiral basis, with ^kx = 
Xck — £f the quasiparticle energy rendered from the 
Fermi level. The off-diagonal elements of the self-energy 
matrix which couples the conduction band and 

valence band, are found to be zero after performing 
the angular integration (this is true also for the pure 
Coulomb interaction case [1|). For TO phonons, Eq. ([6|) 
remains the same except the angular factor becomes 
l=Fcos((/)fe_(-q— 20q) l±cos{(j)k+q — 2(j)q). The Matsubara 



lmtk+{u;) = -^gl (7) 

[^intra(fc, W, UJo)0{uj - UJo)0{uj - UJQ + Ep) 
+-Fintra(fc, t^, -UJq)0{~UJ - UJo)0{uj + UJQ + E p) 
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where the terms Fintra(^, ^^o) = {uj — loq + Ep + 
Ek) (w - Wo + £f + Cfc - |w - Wo + - efcl) /cfc corre- 
sponds to the contribution from quasiparticles making 
an intraband transition in the vicinity of the Fermi 
level by emitting a phonon, and i^intcrlfc, t^, — ^0) = 
- |w + Wo + eF + Cfcl (|t^ + t^o + + ffcl + + t^o + ef - fife) /Cfc 
to quasiparticles making an interband transition from 
the conduction to the valence band by emitting a 
phonon, g1^ = g^AjY^v^ = 2 x 10^^ {A being the area 
of the sample) is the dimensionless phonon-mediated e-e 
coupling constant. Direct calculation of the real part 
of the self-energy from Eq. © gives a logarithmically 
divergent result, which has been noted as well for the 
self-energy calculation with only screened Coulomb 
interaction in the literature This is because the 

linear band structure of graphene in the vicinity of the 
K point is well described by the Dirac Hamiltonian 
Hq, but only up to a cut-off energy scale given by the 
inverse lattice spacing. The real part of the self-energy 
should therefore be regularized by introducing the 
cut-off energy A — hvja'va. the momentum integral in 
Eq. ([5]), whereupon evaluating the integral we obtain 
the following exact analytical expression for ReEfe+(w): 
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In Fig. [T] we show the real and the imaginary parts 
of the self-energy. The gap in ImEfc(w) between the 
phonon energies — wo and wo, characteristic of the op- 
tical phonon emission/absorption process, results from 
the Pauli blocking by electrons located within wo of the 
Fermi level, so that decay by an electron with energy 
w € [— wo,wo] is forbidden because there is no available 
final state to decay to. Beyond the gap, ImEfe(w) be- 
haves linearly as w for large w/wo on both sides of the gap 
due to the linear dependence on energy of the graphene 
density of states. At w < — wq — £fj phonon emission 
can occur through interband transitions, serving as an 
extra decay channel in graphene. The interband contri- 
bution -Fintcr has two dips in |ImI]| below the Dirac point 
(Fig. [1]), one at — wo — £f which comes from phase-space 
restrictions marking the onset of the interband transi- 
tion and the other at — wo — Cfc — £f which originates 
from the angular-dependent electron-phonon interaction 
vertex Eq. ([4]). In regular scmiconductors/metals with 
a parabolic energy dispersion, the imaginary part of the 
self-energy is often calculated under the approximation 
of a constant density of states resulting in a square well- 
shaped profile 6] , which suffices to model the quasiparti- 
cle decay rate quite accurately for these materials. For 
graphene, however, the linear w-dependence of Iml]fe(w) 
away from the gap is a peculiar feature and we believe 
that a correct modeling of the existing experimental data 
for the graphene quasiparticle decay rate ^] should take 
this feature into account. In Fig. [2l we show the cal- 
culated spectral function Afc(w) = — 2ImGfc+(w), where 
Gfc+(w) = l/[w — — Sfe+(w)] is the interacting electron 
Green function with self-energy correction. The spectral 
function comprises two contributions: a delta function 
peak within the gap, and a background beyond the gap. 
As the value of k increases (decreases) , the delta function 
approaches wo (— wo) asymptotically from below (above), 
whereas the background peak at w > (w < 0) becomes 
increasingly prominent. 
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FIG. 2: (Color online) The spectral function Ak{uj) (normal- 
ized to dimensionless unit by multiplying wo) as a function of 
w/wo and k/kp for density n = lO^'^cm"^, the inset shows the 
locus of the position of the quasiparticle delta function inside 
the gap w/wo G [—1, 1]. 



The derivative of the real part of the self-energy 
with respect to w and k determines the effective 
velocity v* , which is given by v*{k)/v ~ (1 + 
B)/{1 - A), where A = 5ReEfc(w)/aw|^^„(fc_fc^) and 

B = {l/v)dReT,k{^^) /dk\^^y(^k~kp)- In metals with a 
parabolic band dispersion, it is a usual practice to ig- 
nore B since in metals ef/^jd ^ 10^ {ojd is the Debye 
energy) and B ^ Yj/ef A ^ S/wd. In graphene, 
however, ef/ojo ~ 1 and B should not be neglected 
in calculating v*/v. Although the phonon-mediated e- 
e coupling (whose magnitude is given by g'^^ ^ 10^^) 
is in general weaker than the Coulomb coupling (given 
by the interaction parameter 0] '"s 0.7), e-ph interac- 
tion actually contributes more significantly to the effec- 
tive velocity renormalization than Coulomb interaction, 
because the real part of the self-energy due to e-ph in- 
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FIG. 3: (Color online) The renormalized conduction band 
energy spectrum for n = 10^''cm~^. The colored intensity plot 
shows the magnitude of the spectral function Ak{u!) whereever 
it is non-zero, while the dashed white line shows the region 
within the amount of phonon energy ujo = 0.196eV from the 
Fermi level where Ak{uj) = 0. The thin dot-dashed straight 
line shows the bare unrenormalized spectrum. 
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FIG. 4: The effective electron-phonon "coupling" parameter 
Acfi = v/v*{kF) — 1 as a function of doping n. 

teraction exhibit sharp changes near the phonon ener- 
gies which are not present in the case of Coulomb in- 
teraction. This results in a larger value of the energy 
derivative of ReE near the phonon energy and therefore 
larger value of v* /v. In Eq. ([8]), logarithmic singu- 
larities occur in ReSfc(a;) at ±wo where Im5]fe(cj) goes 
through a finite step jump, yielding also logarithmic sin- 
gularities in the derivatives of ReSfc(a;) with respect to 
uj and k. Fig. [3] shows the calculated renormalized en- 
ergy spectrum for electron densities n ~ lO^^cm^^, the 
sharp kink shows the logarithmic singularity at w = wq- 
Such a kink only occurs in the conduction band if the 



phonon energy is within the Fermi sea ujq < sp; if the 
phonon energy lies outside of the Fermi sea wq > £f the 
conduction band will be smooth and the logarithmic sin- 
gularity could occur in the valence band. Although our 
theory predicts the correct position of the kink (Fig. [3]) 
in agreement with the experiment 0, the observed en- 
ergy spectrum differs from ours in that the sharpness of 
the kink in the experiment is greatly reduced. First, we 
note that our zero-temperature theory works well in the 
temperature regime of the experiment T < 30K which is 
much smaller than the Fermi temperature Tp ^ 4300K 
at n ~ lO^'^cm^^, hence smearing of the kink due to fi- 
nite temperature effect should be minor. We believe that 
the kink in the experimental spectrum is suppressed due 
to the combined effect of disorder and screening - in par- 
ticular, disorder effects are considerable in the currently 
existing graphene samples. The combination of these ef- 
fects will remove the logarithmic singularity leading to 
a much smoother kink in the energy spectrum, bring- 
ing the calculated spectrum in closer agreement with the 
observed spectrum. Fig. |4] shows the effective electron- 
phonon "coupling" parameter Aeff = v/v*{kp) — 1 as a 
function of electron density n. Our results agree in order 
of magnitude with the extracted value of AoS from the ex- 
periment (Fig. 3c in the second reference of Ref. [il]), and 
the band velocity is shown to be reduced by a percentage 
of {v - v*)/v - 4 - 13% from n = 10" ~ IQi-^cm-^. 

In reality, direct e-e Coulomb interaction, V^^, also 
contributes to the velocity and band renormalization Q 
without, however, contributing any kink structure to the 
energy dispersion. Since both V^J^ and V^^ are weak, 
the leading-order theory including both interactions Q 
would be additive, leading to a Acff = X^g + X^g, where 
Ag^ is shown in Fig. 4 and AJjjj is considered in Ref. [1]. 
We therefore conclude that the graphene velocity renor- 
malization, taking into account both e-ph and Coulomb 
interaction effects, is of the order of 10 — 20% aside 
from the phonon-induccd kink at the phonon energy (i.e. 
Fig. 3). 

In summary, we have developed a theory for the 
phonon-mediated e-e interaction in graphene. We calcu- 
lated the electron self-energy, the spectral function and 
the band velocity renormalization due to the e-ph interac- 
tion. Our results show that the e-ph coupling has a large 
effect on the band structure renormalization, exhibiting 
a kink at ^ 200meV below the Fermi surface as observed 
in the experiment and a reduction of the band velocity 
by ^ 10% — 20% at the experimental doping level. 
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Note added: During the preparation of our manuscript, 
two preprints larXiv:0707.1467t ^l by Calandra et al. and 
arXiv:0707.1666v2 by Park et al., appeared dealing with 
e-ph interaction effects in graphene. 
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